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Abstract:  Using  molecular  dynamics  computer  simulation  technique  we  studied  the 
effect  of  the  ion-electrode  contact  on  the  energetics  of  electron  transfer  reactions.  We 
observed  that  when  the  solvation  shell  of  the  ion  undergoes  large  changes  due  to  the 
ion  transfer  from  the  bulk  to  the  interface,  the  energetics  of  the  reaction  at  the 
interface  differs  substantially  from  the  energetics  of  the  reaction  in  bulk. 
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Introduction. 

*  14  •  . 

Computer  simulations  play  an  important  role'in  testing  the  validity  of  the  ideas 
and  assumptions  made  in  theoretical  description  of  the  observed  phenomena.  For 
example,  in  the  recent  computer  simulation  for  aqueous  ferrous-ferric  electron 
transfer  it  was  observed  that  the  Marcus  relation  [1,2]  for  intersecting  parabolic 
diabatic  free  energy  surfaces  is  quantitatively  accurate  [3]. 

While  many  electron  transfer  processes  occur  in  the  bulk  phase,  some  of  the 
important  electron  transfer  processes  occur  at  the  interface.  Therefore,  it  is  very 
interesting  to  study  how  the  energetics  and  dynamics  of  electron  transfer  process  at 
the  interface  is  different  from  the  process  that  takes  place  in  the  bulk.  Very  recently 
two  simulations  performed  by  Straus  and  Voth  [4]  and  by  Rose  and  Benjamin  [5] 
addressed  this  issue.  Both  simulations  considered  the  electron  transfer  from  the 
surface  of  Pt  metal  to  the  ion  situated  close  to  the  surface.  In  both  simulations  the 
electron  transfer  occurred  when  the  ion  was  situated  ~0.5  nm  away  from  the  Pt 
surface,  meaning  that  the  ion  in  the  simulations  was  still  fully  solvated  by  water,  even 
while  approaching  the  metal  surface.  This  may  correctly  describe  the  situation  for  the 
Fe+3  ion,  but  it  is  not  true  for  an  ion  such  as  I\  Spohr  [6]  and  Perera  and  Berkowitz 
[7]  calculated  the  free  energy  profiles  for  the  I'  ion  solvated  in  water  and  approaching 
the  Pt(100)  surface,  and  obtained  that  most  probable  position  of  this  ion  at  the 
interface  is  when  it  is  in  a  contact  with  the  surface  of  the  metal.  Because  of  the 
strong  interaction  between  Li+  ion  and  Pt  metal  surface,  the  calculation  of  the  free 
energy  profile  for  this  ion  shows  that  its  most  probable  position  is  also  in  a  contact 
with  the  surface  [7].  That  means  that  the  solvation  shells  of  Li+  and  I'  ions  undergo 
a  restructuring  when  these  ions  are  transferred  from  the  bulk  to  the  interface.  The 
influence  of  this  restructuring  on  the  energetics  of  the  electron  transfer  process  from 
and  to  the  ions  is  studied  in  this  paper. 


Molecular  Dynamics  Simulations. 

We  consider  here  four  molecular  dynamics  simulations,  each  performed  in  a 
box,  which  is  1.96*1.96*4.22  nm3  in  dimension  and  which  contains  512  water  molecules 
and  one  ion.  Previous  work  shows  that  a  system  with  these  dimensions  produces  a 
correct  water  density  in  the  center  of  the  box.  The  system  is  periodical  in  x  and  y 
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dimensions  and  is  restricted  by  Pt(lOO)  surface  in  z  dimension.  It  was  also  found  [7] 
that  I*  and  Li*  ions,  when  adsorbed  on  the  surface  of  the  metal,  are  about  0.21  nm 
away  from  plane  of  the  surface.  The  effect  of  the  surface  is  completely  diminished 
when  the  distance  from  the  surface  is  above  1.0  nm.  In  each  of  the  four  simulations 
the  ion  is  kept  at  the  fixed  position.  In  the  first  simulation  (simulation  A)  the  I’  ion 
is  fixed  at  the  interstitial  position  0.21  nm  away  from  Pt  surface.  In  simulation  B  the 
Li*  ion  is  fixed  at  the  same  position.  In  the  simulation  C  the  I'  ion  is  fixed  in  the  bulk 
water  phase  (1.61  nm  away  from  the  surface)  and  in  simulation  D  a  Li*  ion  is  fixed  at 
the  same  distance  from  the  surface. 

The  SPC/E  model  [8]  was  implemented  to  represent  the  water-water 
interaction.  The  ion-water  interaction  was  represented  by  a  Lennard-Jones  plus 
Coulomb  terms,  i.e.  we  used  the  following  form  for  this  interaction: 

UM  -  {At  rn  -  C/r*)+Y'M/rli  (1) 

j 

In  equation  (1)  r  is  the  distance  between  the  ion  center  and  the  oxygen  site  on  the 
water,  q,  the  charge  of  the  ion,  ^  is  the  charge  on  the  site  j  of  water  molecule,  ry  is  the 
distance  between  the  ion  center  and  the  site  j  of  the  water.  The  values  of  the 
parameters  (A  and  C)  for  the  ion-water  Lennard-Jones  interaction  are 
A  *  18927 A12kcal/mol,  C  =  313.2  A6kcal/mol  for  the  Li*  ion  and 
A  *  8993393A12kcal/mol,  C  *  2901.6  A*kcal/mol  for  the  I‘  ion  (In  reference  7,  the 
parameters  A  and  C  for  Li* -water  interaction  were  reported  as  A  *  18927A12kJ/mol, 
C  =  313.2  AakJ/mol,  instead  of  A  *  18927 A12kcal/mol,  C  =  313.2  A6kcal/mol).  To  save 
a  significant  amount  of  computer  time  we  represented  the  water-metal  interaction  in 
the  same  form  as  that  used  previously  [9-11],  This  form  allowed  us  to  concentrate  on 
the  motion  of  water  molecules  only,  replacing  all  the  Pt  atoms  by  an  effective  field. 
The  symmetry  of  the  surface  and  its  corrugation  were  also  taken  into  account  by  the 
effective  potential. 

The  equations  of  motion  were  solved  using  Verlet  algorithm  [12]  and  the  water 
molecules  were  constrained  to  their  geometries  by  implementing  the  Shake  procedure 
[13].  The  timestep  used  in  all  the  simulations  was  2.5  fs.  We  used  a  spherical  cut-off 
scheme,  with  all  interactions  truncated  at  0.98  nm  (a  smoothing  of  the  potential  over 
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the  last  0.05  nm  was  applied).  Corrections  for  the  long  range  interactions  were 
neglected  beyond  the  cut-off  distance.  A  temperature  of  -300  K  was  maintained  by 
the  occasional  rescaling  of  the  velocities. 

Free  Energy  Functions. 

Free  energy  functions  for  the  reactant  and  the  product  are  defined  as: 

Agr(x)  -  -kT  In  PM)  (2) 

Ag„ (X)  -  - kT  In  PM^Gpr  (3) 

where  k  is  the  Boltzmann  constant,  T  is  the  temperature  of  the  system  and  Pr(x)  and 
Pp(x)  are  the  probability  distribution  functions  for  the  reactant  and  product 
respectively,  x  is  a  pre-defined  reaction  coordinate  and  A Gpr  is  the  reaction  free  energy. 
It  was  demonstrated  that  [14,15] 

Agp(x)  -  -x+Ag,(x)  (4) 

and  therefore  we  only  need  to  have  one  smooth  curve  of  either  &gr(x)  or  Agp(x)  to 
obtain  both  curves.  Following  the  well  accepted  convention  [15],  the  reaction 
coordinate  Ae  is  defined  as  the  difference  in  the  solute-solvent  interaction  energy 
between  the  reactant  and  the  product  systems,  i.e. 

Ae  -  er  -  ep  (5) 

For  any  given  configuration,  this  difference  can  be  obtained  in  a  straight  forward 
manner  by  the  following  steps:  first,  calculate  the  interaction  by  setting  the  solute- 
solvent  potential  to  that  of  the  reactant;  then,  calculate  the  interaction  by  setting  the 
potential  function  to  that  of  the  product;  finally,  calculate  the  difference  in  the 
interaction  energies.  In  this  work  we  assume  that  the  Lennard-Jones  part  of  the 
potential  function  is  unchanged  during  the  reaction.  Hence  the  electrostatic  energy 
due  to  the  solute-solvent  interaction  represents  the  reaction  coordinate.  Having 
defined  the  reaction  coordinate  the  probability  distributions  for  the  fluctuations  and 
therefore  free  energy  curves  can  be  calculated  as  functions  of  this  reaction  coordinate. 
But  the  free  energy  curves  can  be  reliably  calculated  only  for  *  sufficiently  close  to  the 
most  probable  value  of  the  probability  distribution  x*  <  Ae>r  or  x=  <  Ae>p.  To  obtain 
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a  smooth  free  energy  curve  in  the  region  far  from  x*  <  Ae>r  or  x=  <  Ae>p  free  energy 
perturbation  and  umbrella  sampling  techniques  can  be  applied  in  the  simulation  [15]. 
A  parameter  a  is  therefore  introduced  to  characterize  the  charging  state  of  the  solute, 

Le.: 


qa  -  (l-o)qr  +  aqp  (6) 

The  solute  has  the  full  charge  of  the  reactant  when  a=0,  it  has  the  charge  of  the 
product  when  a*  1.  A  molecular  dynamics  simulation  is  carried  out  with  the  certain 
value  of  the  parameter  a,  (0<a<  1)  and  the  free  energy  function  is  obtained  from  the 
expression 


A gr°U)  -  -fcTlnP^+ar+AGa,,  (?) 

Increasing  a  stepwise  and  repeatedly  applying  equation  (7),  the  whole  Agr(x)  can  be 
constructed.  In  our  calculations  we  performed  at  least  25  ps  of  molecular  dynamics  for 
each  a  to  get  a  good  convergence.  To  get  a  good  precision  in  the  evaluation  of  solute- 
solvent  energy  fluctuations  at  the  minima  of  free  energy  curves,  the  simulations  at 
these  regions  were  performed  for  at  least  100  ps.  The  curve  for  the  A gp(x)  can  be 
constructed  with  the  help  of  equation  (4).  The  free  energy  difference  A  Ga0  that  enters 
eq.(7)  is  obtained  using  the  following  equation 

AGao  -  -kTln  [<exp-((HcrH0)AT)>0]  <8> 

which  can  also  be  calculated  in  a  stepwise  manner. 

Results  and  Discussions. 

We  performed  our  calculations  for  the  state  when  the  Pt  surface  is  at  the 
potential  of  zero  charge  (pzc),  therefore  all  the  data  for  the  free  energy  of  the  reaction, 
the  transition  free  energy  and  reorganization  free  energy  are  given  with  respect  to  this 
state.  The  calculated  free  energy  curves  Ag/x)  (the  curves  on  the  left)  and  A gp(x)  (the 
curves  on  the  right)  are  shown  in  Figure  la  (for  I'  at  the  interface),  Figure  lb  (for  I' 
in  bulk),  Figure  2a  (for  Li*  at  the  interface)  and  Figure  2b  (for  Li*  in  bulk).  In  Figures 
la  and  lb  the  curves  on  the  left  are  for  the  iodide  ion  in  its  charged  state,  so  that  the 
reaction  corresponds  to  the  process  of  oxidation  of  the  ion.  From  the  comparison  of 
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these  two  figures  we  observe  that  the  partial  removal  of  solvation  water  results  in  a 
smaller  free  energy  of  the  oxidation  reaction,  as  expected.  The  calculations  show  that 
when  the  iodide  is  on  the  surface  AGpr*  153  kJ/mol  and  the  reorganization  energy  is 
195  kJ/mol,  while  when  the  iodide  is  in  the  bulk  AGpr=237  kJ/mol  and  the 
reorganization  energy  is  280  kJ/mol.  The  relatively  large  change  in  the  energetics  of 
the  electron  transfer  reaction  observed  in  the  case  of  I'  ion  (the  reaction  free  energy 
is  changed  by  55%,  while  the  reorganization  energy  is  changed  by  44%)  is  due  to  the 
large  change  in  the  solvation  shell  around  the  ion,  when  it  is  transferred  from  the  bulk 
to  the  surface.  As  it  was  observed  previously  [7]  the  coordination  number  of  I'  changes 
from  ~7  to  ~4  during  this  transfer. 

What  about  the  energetics  of  the  reaction  Li(aq)  -*  Li*(aq)  +  e?  The  free 
energy  curves  for  this  process  are  depicted  in  Figures  2a  and  2b,  where  the  left  curves 
correspond  to  the  uncharged  state  of  the  ion,  so  that  the  reaction  again  corresponds 
to  the  oxidation  process.  The  calculations  show  that  in  case  of  the  lithium  ion  the  free 
energy  of  the  reaction  and  reorganization  free  energy  are  nearly  the  same  when  the 
reaction  takes  place  at  the  surface  or  in  the  bulk.  Thus  for  the  reaction  on  the  surface 
AGp^-512  kJ/mol,  while  for  the  reaction  in  bulk  AGpr=-548  kJ/mol  (the  change  is 
only  7%).  The  absolute  value  of  the  change  of  reorganization  free  energy  is  rather 
large  (reorganization  free  energy  on  the  surface  is  346  kJ/mol,  in  bulk  it  is  equal  to 
445  kJ/mol),  but  the  relative  value  of  the  change  (29%)  is  still  smaller  than  for  I  case 
(44%).  The  smaller  change  in  the  energetics  of  this  reaction  compared  to  the  reaction 
with  iodide  is  explained  by  the  fact  that  the  solvation  of  the  Li+  ion  is  not  changed 
dramatically,  when  the  ion  is  moved  from  the  bulk  phase  to  the  surface.  Thus  Perera 
and  Berkowitz  find  that  although  the  first  solvation  shell  of  water  is  restructured 
around  the  ion,  the  coordination  number  remains  the  same  (~4)  [7]. 

An  issue  often  discussed  in  the  literature  is  how  good  is  the  quadratic 
approximation  to  the  free  energy  curves.  If  the  true  free  energy  curve  is  quadratic, 
then  the  curve  has  a  constant  curvature  that  can  be  obtained  from  the  bottom  region 
of  the  curve.  The  free  energy  curves  with  the  curvatures  obtained  from  such  fit 
("bottom"  curvatures)  are  represented  in  the  Figures  1  and  2  by  thin  lines.  As  we  can 
see  these  curves  approximate  reasonably  well  the  free  energy  curves  for  I'  obtained 
from  the  simulations.  We  also  observe  that  the  "bottom"  curvatures  of  the  reactant 


7 


curves  is  smaller  than  the  "bottom*  curvatures  of  the  product  curves,  i.e.  that  the 
solute-solvent  energy  fluctuations  which  determine  the  magnitude  of  the  curvature 
at  the  bottom  of  free  energy  curves  are  larger  when  the  dynamics  is  performed  with 
the  charged  state  of  the  ion.  This  is  contrary  to  the  behavior  observed  in  the  previous 
simulations  [5].  But  in  the  previous  simulations  the  ion  was  positively  charged,  while 
in  the  present  simulation  the  iodide  ion  is  negatively  charged.  If  the  ion  is  negatively 
charged,  the  neighbouring  water  molecule  places  its  hydrogen  closest  to  the  ion. 
When  the  dynamics  is  performed  on  the  product  surface  the  ion  is  neutral  and 
therefore  the  water  hydrogens  are  further  away  from  the  ion.  Due  to  the  proximity 
of  the  hydrogens  to  the  charged  ion  the  amplitude  of  the  fluctuations  in  ion-water 
energy  is  larger  when  the  dynamics  is  performed  in  the  charged  ion  state. 

As  Figures  2a  and  2b  show  the  quadratic  curves  (with  curvatures  obtained  by 
a  fit  to  the  bottom  of  the  true  free  energies)  do  not  approximate  well  the  curves  from 
the  calculations.  Also  the  "bottom*  curvatures  for  the  reactant  curve  and  the  product 
curves  are  different,  and  again  the  "bottom"  curvature  of  the  reactant  curve  is  smaller 
than  the  "bottom"  curvature  of  the  product  curve.  But  this  time  the  reactant  curve 
is  obtained  from  the  dynamics  with  the  uncharged  ion,  while  the  dynamics  with  the 
charged  ion  is  used  to  get  the  product  curve.  This  time  the  hydrogen  atoms  of  water 
are  closest  to  ions  when  the  dynamics  is  done  with  the  ion  in  its  uncharged  state.  We 
observe  that  the  ion-water  energy  fluctuations  now  have  the  same  character  as  in  the 
previous  simulations  [4,5],  which  were  also  performed  for  the  positively  charged  ions. 

As  was  pointed  out  by  Tachiya  [14],  when  the  "bottom"  curvatures  of  free 
energy  curves  are  different,  the  condition  imposed  by  eq.  (4)  is  violated.  To  satisfy  this 
condition  we  can  fit  one  of  the  free  energy  curves  (for  reactant  or  product)  and  obtain 
the  second  free  energy  curve  from  the  application  of  eq.  (4),  similarly  to  the  procedure 
performed  by  Straus  and  Voth  [4].  The  curves  obtained  in  this  way  are  also  shown 
(broken  lines)  in  Figures  1  and  2.  In  all  cases  we  obtained  the  broken  curves  by 
performing  the  best  fit  to  the  upper  free  energy  curve,  while  the  lower  curve  was 
obtained  from  the  upper  by  using  equation  (4).  As  we  can  see  from  Figures  1  and  2 
this  fit  with  the  "global"  curvature  shows  that  although  we  do  observe  deviations  from 
quadratic  behavior  of  free  energies,  nevertheless  large  regions  of  free  energy  curves 
do  look  parabolic.  And  this  is  correct  for  both  anion  and  cation.  According  to  Marcus 


for  any  intersecting  parabolic  free  energy  curves  of  equal  curvature  there  is  a  universal 
relationship  between  the  activation  free  energy  (AG*)and  the  reaction  free  energy 
(AG°),  which  has  the  form: 

AG  *-  (AG°  +\)2  /4X  (9) 

In  eq.  (9)  X  is  the  reorganization  free  energy.  The  test  of  equation  (9)  for  calculated 
free  energies  is  shown  in  Figure  3.  As  we  can  see  the  calculations  display  deviations 
from  Marcus  curve,  mostly  in  the  inverted  region.  The  closest  to  Marcus  curve 
behavior  is  observed  for  the  reaction  I'(aq)  -*  I(aq)  +  e  in  bulk,  while  the  largest 
deviation  from  Marcus  curve  is  observed  for  the  reaction  Li(aq)  -*  Li+(aq)  +  e  at  the 
interface. 

Figures  4  and  5  show  how  the  reaction  activation  free  energy  depends  on  the 
overvoltage.  We  follow  the  standard  convention  and  assume  that  the  overvoltage  rj  is 
zero  when  the  reaction  free  energy  A is  zero.  Figure  4  presents  plots  for  the 
oxidation  reaction  for  I'  ion  when  it  occurs  on  the  surface  and  in  bulk.  As  we  can  see 
from  this  figure  the  reaction  energetics  on  the  surface  and  in  bulk  is  substantially 
different.  Figure  5  presents  same  type  of  plots  for  the  oxidation  reaction  for  Li.  This 
time  we  do  not  see  a  large  change  in  the  energetics  of  the  reaction  in  agreement  with 
what  was  already  mentioned  above. 

Finally  we  want  to  mention  here  that  the  present  calculations  should  be 
considered  as  model  calculations  rather  than  quantitative,  since  the  quantitative  values 
will  depend  strongly  on  the  potentials  used  in  the  simulations  and  on  the  treatment 
of  long  range  forces  and  the  treatment  of  water  and  ion  polarization.  In  this  paper  we 
used  the  same  ion-surface  interaction  potential  as  the  one  presented  by  Seitz-Beywl 
et  al  in  reference  [16],  where  they  investigated  the  structure  and  dynamics  of  LiF 
aqueous  solution  next  to  Pt(100)  surface.  The  same  potential  was  later  used  by  Spohr 
[6]  and  by  Perera  and  Berkowitz  [7]  to  get  the  potentials  of  mean  force  for  I' 
approaching  Pt(100)  and  Li4  approaching  Pt(100)  surface.  This  ion-surface  potential 
neglects  to  consider  that  partial  charge  transfer  between  the  hydrated  ion  and 
platinum  surface  may  occur.  According  to  estimates  performed  recently  by 
Nazmutdinov  and  Spohr  [17],  partial  charge  transfer  can  be  rather  substantial  in  the 
case  of  I*  ion.  The  ion-surface  potential  that  we  use  also  neglects  the  surface 
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polarization  due  to  the  ion.  The  polarization  of  the  surface  due  to  the  charges  is  often 
treated  using  image  charge  method,  but  this  method  fails  to  describe  appropriately  the 
interaction  of  water  molecules  close  to  the  interface  with  the  surface  of  a  metal  such 
as  Pt  [18].  To  consider  the  effect  of  surface  polarization  due  to  the  ion  on  the 
dynamics  of  water  one  can  use  a  mixed  approach,  where  the  image  charge  method  is 
used  to  describe  ion/metal  interaction,  while  potentials  like  the  one  described  in 
reference  [9]  are  used  to  describe  water/metal  interaction.  Unfortunately,  it  is  not 
known  if  such  an  approach  is  correct.  More  studies  to  test  this  approach  and  to 
investigate  the  effect  of  surface  polarizability  on  the  energetics  and  kinetics  of  electron 
transfer  at  the  interface  are  needed.  Surface  polarizability  may  play  an  important  role 
in  the  electron  transfer  energetics,  but  we  believe  that  the  mqjor  change  in  the 
energetics  (When  such  a  change  exists)  is  due  to  the  large  restructuring  of  the  first 
solvation  shell  when  the  ion  is  adsorbed  on  the  surface,  like  in  the  case  observed  for 
the  I'  ion  in  the  present  simulations. 
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Figure  1.  Free  energy  curves  for  the  reaction,  I'  -*  I  +  e,  (a)  at  the  platinum- 
water  interface;  (b)in  the  bulk  water.  The  reactant  curves,  a  gr(x),  are  on  the  left.  The 
product  curves,  a  gp(x),  are  on  the  right.  The  thick  solid  lines  are  obtained  from  the 
MD  simulations  by  using  Equation  (2),  (3),  and  (7).  The  dashed  lines  are  parabolic 
curves,  k(x  -  b)2  +  c,  fit  to  the  thick  solid  lines.  The  parameters  are  k  =  1.32  x  10"3 
mol/kJ,  b  *46  kJ /mol,  and  c = 154  kJ/mol  for  the  reaction  at  the  interface;  and  k  =  1.00 
x  10*3  mol/kJ,  6=43  kJ/mol,  and  c= 237  kJ/mol  for  the  bulk  reaction.  The  thin  solid 
lines  are  parabolic  curves  with  the  curvature  obtained  from  the  fit  to  the  bottom  of  the 
thick  solid  curves. 

Figure  2.  Free  energy  curves  for  the  reaction,  Li  -*  Li+  +  e,  (a)  at  the 
platinum-water  interface;  (b)  in  the  bulk  water.  Symbols  have  the  same  meanings  as 
in  Figure  1.  The  parameters  are:  k-  3.00  x  10"4  mol/kJ,  b  =  6.12  kJ/mol,  c  =  512 
kJ/mol  for  the  reaction  at  the  interface;  and  k  =  4.66  x  10"4  mol/kJ,  6  =  27.5  kJ/mol, 
and  c  =  548.0  kJ/mol  for  the  reaction  in  the  bulk. 

Figure  3.  The  relationship  between  the  activation  free  energy  -AG*/X  and  the 
reaction  free  energy  -AG°/X.  The  Marcus  relationship  (thick  line)  is  also  plotted  for 
comparison. 

Figure  4.  Activation  free  energy  vs.  the  overvoltage  for  the  reaction,  I*  -» I 
+  e,  in  the  bulk  water  (solid  line)  and  at  the  Pt-water  interface  (dashed  line). 

Figure  5.  Activation  free  energy  vs.  the  overvoltage  ij  for  the  reaction,  Li-*Li+ 
+  e,  in  the  bulk  water  (solid  line)  and  at  the  Pt-water  interface  (dashed  line). 
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